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Since  the  research  project  consists  of  two  parts,  this  final  report  is 
also  divided  into  the  corresponding  two  parts.  The  first  part  summarizes  our 
numerical  results  on  the  flow  of  a  dust-gas  suspension  through  an  oblique 
shock  wave,  and  the  second  part  deals  with  the  flow  of  a  dust-gas 
suspension  having  various  distributions  of  the  properties  of  the  dust 
particles  through  a  normal  shock  wave. 


Abslracl 

The  equations  governing  the  flow  field  which  is  developed  when  a 
supersonic  dustg-gas  suspension  passes  through  a  straight  oblique  shock 
wave  were  formulated.  A  computer  code  for  solving  the  governing  equations 
was  developed  and  used  to  obtain  the  solution  for  a  variety  of  different 
initial  conditions. 

In  addition,  the  dependence  of  the  post-shock  suspension  properties 
on  the  various  physical  properties  of  the  dust  particles,  (namely  the 
diameter  of  the  dust  particles,  their  specific  heat  capacity,  their  material 
density  and  the  loading  ratio  of  the  dust  in  the  suspension)  was 
investigated. 

Nomenclature 

A  -  envelope  area  of  the  solid  particle  =  n02 
C  -  specific  heat  capacity  of  the  solid  particle 
Co  -  drag  coefficient 

Cp  -  specific  heat  capacity  at  constant  pressure  of  the  gaseous  phase 
Cy  -  specific  heat  capacity  at  constant  volume  of  the  gaseo' s  phase 
D  -  diameter  of  the  solid  particle 

Fp  -  drag  force 

k  -  thermal  conductivity 

n  -  flow  Mach  number 

mp-  moss  of  a  solid  particle 
Nu  -  Nusseli  number 

n  -  a  co-ordinate  normal  to  the  oblique  shock  wave 

np  -  number  density  of  the  solid  particles 
P  -  suspension  pressure 

Pr  -  Prandtl  number 
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Qh.t.  *  rate  of  heat  per  unit  volume  transferred  from  the  gaseous  to  the  solid 
phase 

R  -  specific  gas  constant  > 

Re  '  Reynolds  number 

5  -  cross  section  of  the  solid  particle  =  nD^M 
T  -  temperature 

u  -  velocity  of  gaseous  phase 
V  -  velocity  of  the  solid  phase 

y  -  specific  heat  capacities  ratio  s  Cp/Cv 
^  ‘  angle  of  incidence 

6  -  flow  deflection  angle 
p  -  dynamic  viscosity 

p  '  spatial  density 
a  '  material  density 

T)  -  loading  ratio  of  the  solid  phase  in  the  suspension  :  Ppq/pqq 

Subscripts 

g  -  gaseous  phase 

n  ‘  normal  component 

p  -  solid  phase 

s  -  tangential  component 

X  -  horizontal  co-ordinate 

y  -  vertical  co-ordinate 

0  -  flow  state  ahead  of  the  shock  wave 

t  -  flow  stale  behind  the  shock  wave 

Introduction 

The  interest  in  the  gas-dynamic  behaviour  of  a  gas-particle 
suspension  grew  in  the  past  three  decades  due  to  its  application  to  many 
engineering  problems.  Some  typical  examples  are:  metallized  propellents  of 
rockets,  jet-type  dust  collectors  end  blast  waves  in  dusty  atmospheres.  In 
addition,  mixtures  with  gases  heavily  laden  with  particles  occur  frequently 
in  industrial  processes  such  as  plastics  manufacturing,  flour  milling,  coal- 
dust  conveying,  powder  metallurgy  and  powdered-food  processing.  General 
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descriptions  of  such  flows  can  he  found  in  several  books  end  review  papers 
ISoo  ( 1 967),  Marble  ( 1 970)  and  Rudinger  ( 1 973)] 

The  major  differences  between  the  flow  fields  which  are  developed 
behind  a  normal  shock  wave  in  a  dusty*gas  and  a  pure  (dust*free)  gas  are 
illustrated  in  figures  la  and  b  for  the  temperatures  and  the  velocities, 
respectively.  When  a  steady  pure  gas  encounters  a  normal  shock  wave  it 
experiences  a  sharp  (almost  discontinuous)  change  in  its  thermodynamic  and 
kinematic  properties.  This  sudden  change  is  shown  in  figure  I  to  occur 
between  (0)  and  (I).  The  thickness  of  this  disturbance.  If,  is  only  a  few 

mean  free  paths  of  the  gas  atoms  or  molecules  (about  6.6  x  10'^  cm  in 
standard  conditions).  Beyond  (1)  the  gas  properties  remain  constant  (solid 
lines  in  figures  la  and  b)  provided  the  gas  conditions  at  (1)  are  not 
sufficient  to  excite  the  internal  degrees  of  freedom  of  the  gas. 

If,  however,  the  gas  is  laden  with  solid  particles  then  the  suspension 
which  was  originally  at  a  state  of  thermodynamic  and  kinematic 
equilibrium,  ahead  of  the  shock  front,  is  suddenly  changed  into  a  non* 
equilibrium  state,  because  the  solid  particles,  due  to  their  size  compared 
with  If,  do  not  experience  any  noticeable  change  in  their  properties  upon 

moving  from  (0)  to  (I).  Thus,  at  (1)  the  gas  has  a  much  higher  temperature 
than  the  dust,  Tg  >  >  Tp  and  a  much  lower  velocity  u  <  <  v  .  Morgenthaler 
(1962)  indicated  that  this  is  true  even  if  the  particle  diameter  is  as  small 
as  0.1  pm  (for  shock  waves  in  air  at  nearly  standard  conditions,  where  the 
mean  free  path  is  about  0.066  pm).  Therefore,  the  particles  are  not 
influenced  by  the  initial  disturbance,  and  the  gas  properties  at  (I)  can  be 
safely  assumed  to  be  identical  to  those  of  a  pure  gas  with  the  some  initial 
conditions. 

Far  downstream  of  (1),  i.e.,  at  (m)  in  figure  I,  the  gas  and  the  solid 
phases  reach  a  new  state  of  thermodynamic  and  kinematic  equilibrium  via 
momentum  and  energy  exchange.  Theoretically  all  shock  waves  in  dusty 
gases  are  infinitely  thick,  since  equilibrium  is  approached  asymptotically. 
However,  it  is  a  common  practice  to  assign  to  the  shock  wave  an  effective 
thickness  which  is  defined  by  a  requirement  that  the  suspension  properties 
come  close  to  their  equilibrium  downstream  values.  It  was  shown  by 
Gottlieb  &  Coskunses  (1985)  that  the  suspension  equilibrium  properties  (at 
infinity)  can  be  calculated  from  the  usual  normal  shock  wove  relations. 
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provided  that  the  usual  pure  gas  parameters  y  and  R  (the  specific  heat 
capacities  ratio  and  the  specific  gas  constant)  are  replaced  by  effective 
values  y  ahd  R  which  solely  depend  04  the  initial  conditions  of  the 
suspension. 

Between  (1)  and  (m)  the  solid  particles  are  not  in  equilibrium  with 
the  gas.  The  flow  region  between  (I)  and  (m)  is  known  as  the  relaxation 
zone«  for  it  is  analogous  to  the  relaxation  zone  in  pure  gases  where  the 
internal  degrees  of  freedom  are  excited.  The  extent  of  the  relaxation  zone 
strongly  depends  on  the  momentum  and  heat  transfer  mechanisms  which 
enable  the  solid  and  the  gaseous  phases  to  reach  a  new  equilibrium  state. 
The  analysis  of  the  relaxation  zone  was  studied  by  many  investigators.  The 
pioneering  works  of  Carrier  (1958),  Kriebel  (1964)  and  Rudinger  (1964) 
verified  the  exist^^  .e  of  this  relaxation  zone  and  identified  the  parameters 
affecting  it,  namely;  the  solid  particle  diameter,  D,  its  specific  heat 
capacity,  C,  its  material  density,  a,  and  the  loading  ratio,  q.  Igra  &  6en*Dor 
(1980)  compared  various  correlations  for  the  drag  coefficient,  Cq,  and  the 

heat  transfer  coefficient,  Nu,  and  pointed  out  their  effect  on  the  extent  of 
the  relaxation  zone.  In  addition  they  studied  the  role  of  thermal  radiation 
heat  transfer  between  the  two  phases  and  showed  that  it  can  be  neglected 
when  the  incident  shock  waves  tiach  number  is  smaller  than  five. 

in  all  the  above  mentioned  works,  as  well  as  in  many  others,  the 
gaseous  phase  was  assumed  to  behove  as  0  perfect  gas.  This  assumption 
was  relaxed  by  Ben-Dor&  Igra  (1982)  and  Igra  &  Ben*Dor  (1984)  who  solved 
the  flow  field  while  accounting  for  real  gas  effects.  Dissociating  nitrogen 
was  the  gaseous  phase  in  the  tetter  work  end  ionizing  argon  in  the  former. 

The  assumption  that  the  solid  particles  are  inert,  which  was  also 
adopted  in  most  of  the  published  studies,  was  relaxed  by  Elpehn,  6en-0or  & 
Igra  (1986)  who  solved  the  flow  field  of  an  oxygen*cerbon  suspension 
passing  through  a  normal  shock  wave,  behind  which  the  temperature  of  the 
carbon  particles  reached  their  ignition  temperature  and  burned  out. 

The  assumption  of  uniform  solid  particles  was  relaxed  by  Elata,  Ben- 
Oor&  Igra  (1988)  who  solved  the  cose  of  size-distributed  solid  particles. 

In  all  the  above  mentioned  solutions  the  flow  field  was  one¬ 
dimensional  end  steady.  The  aim  of  the  present  study  was  to  solve  the  case 
of  a  two-dimensional  steady  flow.  This  is  the  case  when  the  shock  wave  is 
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oblique.  There  are  many  incidences  where  the  shock  wave  is  oblique.  For 
example,  one  can  mention  the  shock  wave  generated  by  a  supersonic  vehicle, 
the  shock  wave  which  is  developed  at  the  entrance  nozzle  of  a  rocket  engine 
and  the  reflected  shock  waves  which  arise  when  an  explosion  generated 
blast  wave  interacts  with  man*made  structures.  In  all  these  cases  the 
shock  wave  is  oblique  and,  hence,  unlike  the  previously  mentioned  cases  the 
resulted  flow  field  behind  the  shock  wave  is  two-dimensional. 

The  aim  of  the  present  study,  therefore,  is  to  solve  the  flow  field 
which  is  developed  when  a  supersonic  dusty  gas  suspension  passes  through  a 
straight  oblique  shock  wave. 

Figure  2  illustrates  schematically  the  problem  to  be  solved.  A  dust 
gas  suspension  which  is  in  a  thermal  and  kinematic  equilibrium  encounters 
an  oblique  shock  wove.  The  angle  of  incidence  is  i  (sometimes  known  as  the 
wave  angle).  As  mentioned  previously,  upon  the  passage  of  the  suspension 
through  the  shock  wave,  the  properties  of  the  gaseous  phase  assume  a  new 
state,  known  as  the  frozen  state,  immediately  behind  the  shock  front,  while 
the  solid  phase  passes  through  the  shock  wave  unaffected.  Due  to  the  fact 
that  the  shock  wave  is  oblique,  the  streamline  of  the  gaseous  phase  is 
deflected  by  an  angle,  6g,  while  the  trajectory  of  dust  particles  remains 
unchanged.  Thus,  in  addition  to  the  differences  in  the  properties  of  the  two 
phases  which  were  mentioned  previously,  here  there  is  also  a  difference  in 
the  direction  of  propagation  of  the  two  phases.  The  two  phases  which  are  no 
longer  in  equilibrium,  start  exchanging  momentum  and  energy  until  they 
finally  reach  a  new  thermal  and  kinematic  equilibrium.  Here,  kinematic 
equilibrium  means  that,  in  addition  to  equal  velocities,  the  two  phases  also 
reach  the  same  direction.  The  flow  region  inside  which  the  two  phases  are 
in  non-equilibrium,  i.e.,  the  relaxation  zone,  is  also  shown  in  figure  2. 

Theoretical  Background 

Assumptions 

The  assumptions  upon  which  the  present  model  is  based  and  their 
implications  ere  listed  in  the  following: 

I)  The  gaseous  phase  behaves  as  an  ideal  gas.  Thus,  the  equation  of  stole 
of  the  gas  Is  P  =  p^RTg.  Note  that  It  Is  not  assumed  here  that  the  gas  is 
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calohcatly  ideal.  Alternatively,  the  dependence  of  both  Cp  and  Cy  on 

the  gas  temperature  is  accounted  for.  This  has  not  been  done  in 
previous  studies  where  both  Cp  and  Cy  were  assumed  to  be  constant. 

2}  All  solid  particles  are  rigid,  inert  small  identical  spheres  uniformly 
distributed  in  the  gaseous  phase.  Thus  there  is  no  heat  addition  or 
reduction  due  to  chemical  processes  between  the  solid  and  the  gaseous 
phases.  Furthermore,  Re  and  Nu  are  based  on  the  particle  diameter,  0. 

3)  The  volume  of  the  solid  phase  in  the  suspension  can  be  neglected.  Thus 
the  momentum  end  energy  exchange  among  the  solid  particles  can  be 
ignored. 

4)  Aside  from  momentum  and  energy  interactions  between  the  gaseous  and 
the  solid  phases,  the  gaseous  phase  is  considered  to  be  a  perfect  flow, 
i.e.,  the  dynamic  viscosity,  |i^,  and  the  thermal  conductivity,  are 
zero.  This  also  implies  that  neither  kinematic  nor  thermal  boundary 
layers  develop  around  the  solid  particles. 

5)  The  dynamic  viscosity,  the  thermal  conductivity,  kg,  and  the 
specific  heat  capacity  at  constant  pressure,  Cp,  of  the  gaseous  phase 
depend  solely  on  its  temperature,  Tg. 

6)  The  solid  particles  are  too  large  to  experience  any  change  in  their 
thermodynamic  and  dynamic  properties  upon  their  passage  through  the 
shock  front.  In  addition  they  are  also  large  enough  not  to  experience 
Brownian  motion.  Thus,  the  partial  pressure  of  the  solid  phases  can  be 
neglected. 

7)  The  solid  particles  are  small  enough  to  satisfy  the  condition  Bj  <  0.1, 
where  Bj  is  the  Biot  number,  Bj  =  hr/kp  (h  is  the  coefficient  of 
connection  heat  transfer,  r  is  the  radius  of  the  solid  particle,  and  kp  is 
its  thermal  conductivity).  Thus  the  temperature  within  the  solid 
particles  can  be  assumed  to  be  uniform. 

B)  The  weight  of  the  solid  particles  and  the  buoyancy  forces  experienced 
by  them  are  negligibly  small  in  comparison  with  the  drag  forces  acting 
on  them. 

9)  The  specific  heat  capacity,  C,  of  the  solid  particles  is  constant. 

10)  Ahead  of  the  normal  shock  wave  the  suspension  is  at  a  state  of 
thermodynamic  and  kinematic  equilibrium,  i.e.,  uq  =  vq  and  T^q  =  Tp^. 

where  u  and  v  are  the  velocities  of  the  gas  and  the  solid  particles,  and 
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Tg  and  Tp  are  the  temperatures  of  the  gaseous  and  solid  phases, 
respectively. 

1 1}  Based«n  the  density  ratio  of  the  two  phases  1/2500.  the  virtual  mass 
which  depends  on  this  ratio  is  neglected. 

In  addition  to  the  above  listed  assumptions  it  is  assumed  that  the  flow  iielo 
under  consideration  is  two-dimensional  and  steady. 

if  the  entire  problem  is  analyzed  in  the  (s.n)-plone.  where  the  s-oxis 
is  parallel  to  the  shock  wave  front  and  the  n-axis  is  normal  to  the  shock 
wave  front,  then  the  problem  at  hand  can  be  considered  as  one-dimensional. 
Note  that  the  s-component  of  the  velocities  of  the  two  phases  immediately 
behind  the  shock  front  ore  identical,  i.e..  uij  =  vi).  For  this  reason  the 
relative  velocity  between  the  two  phases  in  the  s-direction  is  zero,  and 
therefore  there  are  no  drag  forces  in  the  s-direction.  As  a  consequence, 
there  are  no  changes  in  the  velocity  components  in  the  s-direction  and  they 
remain  constant  in  the  entire  flow  field.  From  symmetrical  considerations, 
it  is  obvious  that  all  the  properties  of  the  suspension  must  remain  constant 
along  n=constant  lines,  i.e..  along  lines  which  are  parallel  to  the  s-axis. 
Therefore,  the  partial  derivative  of  any  of  the  suspension  properties  with 
respect  to  s  is  zero. 

Governing  Equations 

Based  on  the  foregoing  discussion  and  assumptions,  the  governing 
equations  for  the  problem  at  hand  are: 

-  conservation  of  mass  of  the  gaseous  phase: 

“^(pgUn)  =  0  (1) 


-  conservation  of  mass  of  the  solid  phase: 

-|;(ppVn)  =  0 


(2) 


a 


-  conservation  of  linear  momentum  of  the  gaseous  phase: 


c 

dn  ■ 


'(3) 


-  conservation  of  linear  momentum  of  the  solid  phase: 


PpVn~~=-Foo  (4) 

-  conservation  of  energy  of  the  gaseous  phase:  - 
d  1  2 

PgUn  (CpTg  ♦  “2  u^)  =  Qh.T.  ♦  FOnVn  (5) 

>  conservation  of  energy  of  the  solid  phase: 
d  1  2 

Pp'^n  ^^P  “2  s  -Qh.t.  -  FpnVn  (6) 

where  pg  and  pp  are  the  spatial  densities  of  the  gaseous  and  solid  phases,  un 
and  Vn  are  the  velocity  components  of  the  gaseous  and  solid  phases  in  the  n- 
direction,  Tg  and  Tp  are  the  temperatures  of  the  gaseous  and  solid  phases,  P 
is  the  pressure  of  the  suspension,  Cp  is  the  specific  heat  capacity  at 
constant  pressure  of  the  gaseous  phase,  C  is  the  specific  heat  capacity  of 
the  solid  particles.  Qh.t  is  the  amount  of  heat  transferred  per  unit  volume 
from  the  gaseous  phase  to  the  solid  phase  and  Fpn  is  the  drag  force  per  unit 
volume  applied  by  the  gaseous  phase  on  the  solid  particles. 

Thus,  if  Cp,  C,  Qh.t  and  Fpn  ore  expressed  in  terms  of  the  dependent 
variables  namely;  pg.  pp,  Un,  Vn,  Tg.  Tp  and  P,  then  the  above  set  of  six 
conservation  equations  contain  seven  unknowns.  The  seventh  equation  which 
is  required  to  make  the  above  set  of  equations  solvable,  is  the  equation  of 
stole  of  the  gaseous  phase;  i.e.. 


P  =  pg  R  Tg 


(7) 
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In  the  following^  the  above  mentioned  complementary  equations  are 
developed. 

The  drag  force  per  unit  volume  can  be  obtained  by  multiplying  the 
drag  force  acting  on  a  single  solid  particle  by  the  number  density  of  the 
solid  particles. 


Fon  =  "2  Pfl  (Vn“Uft)  Ivn-Uftl  S  np 


(8) 


where  Cp  is  the  drag  coefficient,  S  is  the  projected  cross  section  of  the 
solid  spherical  particle  (i.e.,  S  =  n02/4,  where  D  is  the  diameter  of  the  solid 
particles),  and  np  the  number  density  of  the  solid  particles  can  be  calculated 
from 


np  = 


-£l 

mp 


(9) 


where  mp  is  the  mass  of  a  single  solid  particle  (i.e.,  mp  s  onOVd.  where  a  is 
the  material  density  of  the  solid  particles). 

The  drag  coefficient,  Cp,  is  usually  expressed  as  a  function  of  the 
Reynolds  number,  Re,  which  in  turn  depends  on  the  slip  velocity  Ivn-uni.  For 
this  reason,  the  Re  number  is  very  high  immediately  behind  the  shock  wave 
while  its  magnitude  vanishes  towards  the  end  of  the  relaxation  zone,  where 
vn  s  Un,  is  approached.  For  this  reason  two  different  correlations  for  the 
drag  coef  f icient  Cp  ere  used: 

For  Re  <  600 


Cp  =  “~(l  ♦O.I5ReO-687) 

and  for 800  <  Re  <  3x  lOS 

24  0  42 

*  <^  l5ReO-687)^  ,^42500  Re*l  >8 


(lOa) 


(100) 


60th  of  these  correlations  were  taken  from  Clift.  Grace  &  Weber  (1976).  The 
Reynolds  number  Is  calculated  from: 
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Re , 

i  H 


(n) 


where  the  dynamic  viscosity  of  the  gaseous  phase,  is  calculated  t>y  the 
expression  suggested  by  hazor.  Ben'Oor&  Igra  (1985)  which  reads: 

Where  is  the  dynamic  viscosity  of  the  gaseous  phase  at  a  reference 
temperature  Tq^. 

The  heat  transferred  per  unit  volume  from  the  gaseous  phase  to  the 
solid  phase  can  be  obtained  by  multiplying  the  heat  transferred  to  a  single 
particle  by  tho  number  density  of  the  solid  particles,  i.e.. 

Qh.t  =  Ah(Tp“T9)np  (13) 

where  A  is  the  surface  area  of  the  solid  spherical  particle  (i.e..  A  s  nO^)  and 
h  the  coefficient  of  heat  convection  can  be  calculated  from 


h  = 


Nukq 

0 


(14) 


where  is  the  thermal  conductivity  of  the  gaseous  phase  and  Nu  the 
Nusselt  number  is  a  function  of  the  Reynolds  number.  Re.  and  the  Prandtl 
number.  Pr.  which  is  defined  as 

Pr-.-^  (IS) 

k. 

The  correlation  for  Nu  which  was  used  in  the  present  study  is  taken 
from  Drake  (1961).  It  reads; 


Nu  =  2  ♦  0.459  Prl/3  ReO-SS 


n 


As  mentioned  in  assumption  (S)  the  specific  heat  capacity  at  constant 
pressure^  Cp,  and  the  thermal  conductivity^  depend  solely  on  the 
temperature>of  the  gaseous  phase,  Tp.  For  this  reason  the  gaseous  phase  was 
identified  as  air.  For  air  (Holman,  1981,  p.  542)  in  the  range 
300<Tfl<2500K 

kg  =  0.0125  ♦  5.509  X  10-5  Tg  (16) 

and  from  Von  Wylen  &  Sonntag  (1978)  p.  683 

C,  =  39.06  -  512.79  '  ^  ♦  1072.7  ^  -  820.4 


this  condition  is  good  for  300  <  Tg  <  3500  K. 

Numerical  Solution 

Now  we  are  at  the  stage  where  the  governing  equations  are  well 
defined,  i.e.,  the  number  of  independent  variables  is  equal  to  the  number  of 
equations  and  all  the  non*dependent  variables  can  be  expressed  in  terms  of 
the  dependent  variables  and  other  known  parameters. 

As  can  be  seen,  the  governing  equations  of  the  problem  at  hand 
consist  of  seven  non-linear  differential  equations  which  must  be  solved 
simultaneously  in  order  to  obtain  the  spatial  distribution  of  the  flow 
properties. 

Numerous  computer  code  packages,  capable  of  numerically  solving 
non-linear  differential  equations,  are  available,  e.g.,  the  IMSL  (international 
hothematical  end  Statistical  libraries).  This  package  contains  three 
different  computer  codes  for  solving  a  given  set  of  differentia)  equations 
provided  the  initial  conditions  ore  known.  (As  will  be  shown  subsequently, 
the  initial  conditions  in  the  case  at  hand  ere  indeed  known.)  The  most 
accurate  code  out  of  these  three  codes  is  DREBS.  It  is  based  on  the 
extrapolation  method,  and  Is  superior  when  very  high  accuracy  is  required 
and  when  the  calculation  of  the  derivatives  is  not  expensive. 

Fortunately,  ^  ie  relative  simplicity  of  the  governing  equations  of  the 
problem  at  hand,  enables  one  to  rewrite  them  in  a  form  where  the 


derivatives  of  the  dependent  variables  are  isolated,  and  hence  can  be 
calculated  very  quickly  and  cheaply  (computer-wise).  The  rewritten  set  of 
the  governing  equations  con  be  changed  to  the  following  form:  i 


>/  dTg  Td  du\ 

■sr=  p’®  (to  u  "toj 

dpq  pq  dU 

dn  '  ”  u  dn 

dtp  Qht. 
dn  *  cppv 

dv  Fq 
dn  ~  ppv 

dpp  ^  pp  dv 
dn  ”  *  v  dn 


(17) 

(18) 

(19) 

(20) 

(21) 

(22) 

(23) 


The  Initial  Conditions 

As  mentioned  earlier,  the  problem  at  hand  is  such  that  the  initial 
conditions,  i.e..  the  values  of  the  dependent  variables  at  k  s  0  (at  stale  (l)  in 
figure  1 1  can  be  easily  calculated. 

Based  on  the  fact  that  the  thickness  of  the  shock  wave  front  is  orders 
of  magnitude  smaller  than  the  diameter  of  the  solid  particles,  it  is  a 
common  practice  to  assume  that: 

I.  The  gaseous  phase  experiences  the  well  known  ‘frozen’  change  upon  its 
passage  through  the  shock  wove. 
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2.  The  solid  phase  remains  unaffected  as  it  passes  through  the  shock  wave. 

The  above  assumptions  imply  that  while  the  properties  o^  the  solid 
phase  immediately  behind  the  shock  wave  front  ore  identical  to  jihose  ai.sad 
of  it,  the  properties  of  the  gaseous  phase  can  be  calculated  using  the  well 
known  Rankine-Hugoniot  relations  which  relate  the  gas  properties  in  both 
sides  of  the  shock  wave.  More  specifically; 


P9|  =Pso 


(rO  Mq  *2^ 


Pi  =Po 


(24) 


(25) 


"oX- 


TT”o 


iulE 

2(ri)  0 


(26) 


ni  = 


2 

xl 


-  1 

r  I  0 


Uf  =  rii  RTgj 


(27) 


(28) 


where  subscripts  '0*  and  ‘I*  denote  the  flow  states  immediately  ahead  and 
behind  the  shock  wave.  f1  is  the  Mach  number  of  the  gaseous  phase,  i.e., 
tl  s  u/a  (a,  the  local  speed  of  sound,  is  simply  calculated  from  a  =  (y  R  T)i/2). 

Note  that  when  Mo  i.e.,  when  the  shock  wave  degenerates  to  a 
Mach  wave,  then  there  is  no  change  in  the  properties  of  the  gaseous  phase  as 
it  passes  through  the  shock  wave. 

Since  the  model  at  hand  assumes  that  the  gaseous  phase  behaves  as  a 
perfect  gas  the  upper  value  of  ho  which  can  be  used  is  limited,  since  as 
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Mo  ->  M  Tf  M.  In  the  case  of  a  diatomic  gas,  such  as  oxygen  (O2),  nitrogen 
(N2),  etc..  Mo  s  6  is  usually  taken  as  the  upper  limit  for  which  the 
assumption  of  a  perfect  gas  behaviour  is  valid.  For  monatomic  gapes  such  as 
helium  (He),  neon  (Ne),  argon  (Ar),  etc.,  the  upper  limit  of  Mo  is  even  higher. 
Beyond  these  limiting  values,  reel  gas  effects  must  be  accounted  for.  The 
ways  of  treating  such  cases  are  discussed  by  Ben-Dor  and  Igra  (1982)  and 
Igra  and  Ben-Dor  (1984). 

Complementary  Equations 

Once  the  set  of  the  governing  equations  is  solved,  the  velocity 
vectors  of  each  of  the  two  phases  can  be  calculated  from: 

I 

u  =  Uft  n  ♦  u*  s  (29a) 

V  s  Vft  n  ♦  Vj  s  (29b) 

The  absolute  velocities  of  the  two  phase  ore  therefore: 

|u|  =  (u^^u^)>/2  (308) 

n  9 

Iv|  =  (v^tv^)'f2  (30b) 

n  9 

The  deflection  of  the  gaseous  phase  streamline,  Sg,  from  its  original 
direction  can  be  calculated  from 

eg  =  i-tan-l-{7  (31a) 

Similarly,  the  deflection  of  the  dusty  phase  streamline.  Bp,  is 


0p  =  ♦  -  tan"  1 


-1  Jin 

V9 


(31b) 
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Obviously,  when  the  suspension  reaches  a  new  equilibrium  state  at  the  end 

of  the  relaxation  2one  the  flow  directions  of  the  two  phases  become  parallel 
and  hence  s  Bp^.  ^ 

The  velocity  components  of  the  gaseous  and  the  dusty  phases  in  the 
(x,y)'Coordinate  system  shown  in  Figure  I  can  be  calculated  from: 


Ux  =  Un  sin^  ♦  Us  COS4  (32a) 

Uy  s  -Un  cos^  *  Us  sin^  (32b) 

and  similarly  for  the  solid  phase: 

Vx  =  Vn  sin*  ♦  Vs  COS4  (33a) 

Vy  =  -vn  cos*  ♦  Vs  sin*  (33b) 


Numerical  Results 

As  mentioned  earlier,  the  problem  at  hand  can  be  treated  as  one¬ 
dimensional  by  solving  it  in  the  (s,n)-plane.  For  this  reason,  the  flow 
profiles  shown  in  Igra  and  Ben-Dor  (1960),  for  example,  are  all  applicable 
for  the  case  of  oblique  shock  waves,  provided  the  flow  Mach  number  Hq  is 
replaced  by  no  sin*  and  the  x-axis  is  replaced  by  the  n-axis.  Note  that  if  the 
angle  of  incidence  is  set  to  *  s  90^,  then  the  oblique  shock  wave  becomes 
normal  to  the  oncoming  flow,  and  the  general  case  of  an  oblique  shock  wave 
degenerates  to  the  well  known  one-dimensional  normal  shock  wove  case. 

Since  the  change  of  the  flow  properties  in  a  direction  normal  to  the 
oblique  shock  wave  front  as  well  as  the  dependence  of  these  properties  on 
the  physical  properties  of  the  dust  particles  can  be  deduced  from  Igra  and 
Ben-Dor  (I960),  in  the  following  only  results  which  are  unique  to  the  fact 
that  the  shock  wove  is  oblique  are  presented. 

The  Gas  Deflection  -  By 

As  mentioned  earlier,  the  direction  of  the  streamline  of  the  gaseous 
phase  changes  continuously  from  its  'frozen*  direction  to  its  'equilibnum* 
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direction,  in  the  following,  the  dependence  of  the  flow  deflection  ongle,  6^, 
on  the  physical  properties  of  the  dust  is  discussed. 

The  dependence  of  the  flow  deflection,  6q,  on  the  dust  Iqoding  ratio, 
n.  is  illustrated  in  figure  3,  for  Mo  s  3  and  ^  =  30<3.  The  dust  physical 
parameters  are  D  =  lOO^im,  C  s  tOOOJ/CKgmK),  a  =  1500Kgm/m3. 

As  can  be  seen  the  frozen  value  of  6),  i.e.,  6^  at  n  =  0  is  identical  for 
all  the  cases.  However,  the  larger  the  loading  ratio  is,  the  greater  the 
equilibrium  flow  deflection  becomes.  The  dashed  lines  are  the  values 
calculated  by  the  equivalent  gas  concept  presented  in  the  introduction, 
which  should  be  reached  at  the  end  of  the  relaxation  zone. 

It  is  also  evioent  from  figure  3  that  the  larger  the  loading  ratio  is, 
the  shorter  the  relaxation  length  becomes.  While  for  n  =  O.l  it  takes  about 
am  for  the  gaseous  phase  to  reach  the  equilibrium  deflection  angle,  only  3m 
are  required  when  q  =  2. 

The  deflection  angle  at  the  end  of  the  relaxation  zone,  as  a 

function  of  the  dust  loading  ratio,  q,  is  shown  in  a  different  way  in  figure  A 
for  Mq  =  d.  ♦  =  600, 0  -  loOiim,  C  a  1000J/(KgmK)  and  a  s  iSOOKgm/m^. 

The  fact  that  the  equilibrium  deflection  angle  increases  with 
increasing  loading  ratios  is  clearly  seen  in  figure  A.  However,  it  is  evident 
that  as  q  increases  the  rate  of  increase  of  the  equilibrium  deflection  angle 
decreases  end  it  seems  that  there  is  an  upper  limit  on  69^^  as  q  assumes 

very  high  values. 

The  dependence  of  the  deflection  angle,  69,  on  the  diameter  of  the 
solid  particles  is  shown  in  figure  5.  Again  do  =  3  and  ^  =  300.  The  dust 
properties  are  q  =  0.5,  C  s  1000J/(Kgm  K)  and  0  =  ISOOKgm/m^. 

Since  the  loading  ratio,  q,  and  the  specific  heat  capacity  of  the  dust 
particles,  C,  are  the  same  for  the  cases,  the  equilibrium  values  at  the  end  of 
the  relaxation  zone  ere  identical.  However,  it  is  evident  from  figure  5  that 
the  smaller  the  diameter  of  the  solid  particle  is,  the  shorter  the  relaxation 
length  becomes. 

The  dependence  of  the  deflection  angle,  B9,  on  the  specific  heat 
capacity  of  the  solid  particles,  C,  for  hq  =  3  and  ^  =  30^  is  shown  in  figure 
6.  The  properties  of  the  solid  phase  ore  q  :  0.5,  D  =  lOOpm  and 
o  =  l500Kgm/m3  It  is  evident  from  figure  6  that  higher  values  of  the 
specific  heat  capacities  result  in  larger  deflections  at  the  end  of  the 
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relaxation  zone.  The  length  of  the  relaxation  zone  is  not  seen  to  strongly 
depend  on  the  specific  heat  capacity. 

The  dependence  of  the  deflection  angle,  6g,  on  the  material  density  oi 
the  solid  particles  for  Mq  =  3  and  ♦  =  30^  is  shown  in  figure  7.  The  dust 
properties  are  q  =  0.5,  D  =  lOO^ifn,  C  =  1000J/(Kgm  K).  Again,  since  n  and  C 
are  Identical  for  ell  the  cases,  the  equilibrium  values  approached  towards 
the  end  of  the  relaxation  zone  ore  the  same  for  all  the  calculated  cases,  it  is 
clearly  seen  from  figure  7  that  the  smaller  the  matenal  density  is.  the 
shorter  the  relaxation  zone  becomes.  While  it  is  about  4m  for 
o  =  lOOOKgm/m^  it  increases  to  about  lOm  when  a  is  increased  to 
2500Kgm/m5. 

The  gas  and  the  solid  particle  paths  is  shown  in  figure  6  for  Mq  =  3. 

<t>  =  600,  q  =  1.0,  D  =  lOO^im,  C  =  IOOOJ/(Kgm  K)  and  a  =  1500Kgm/m3.  The 
difference  between  the  trajectories  is  clearly  evident.  While  the  gas 
particles  are  deflected  immediately  upon  their  passage  through  the  shock 
wave  (n=0),  the  solid  particles  do  not  change  their  original  direction. 
However,  behind  the  shock  front  the  streamlines  of  the  two  phases  approach 
a  parallel  direction.  The  results  shown  in  figure  8  indicate  that  it  is 
impossible  to  shape  a  wedge  in  a  supersonic  flow  of  a  dusty  gas,  in  such  a 
way  that  it  will  generate  a  straight  oblique  shock  wave.  This  is  due  to  the 
fact  that  a  straight  oblique  shock  wave  in  a  dusty  gas  results  in  two 
different  particle  trajectories,  and  no  wedge  can  satisfy  these  two  profiles. 
However,  if  the  wedge  surface  is  sticky,  i.e.,  if  it  con  be  assumed  that  when 
a  solid  particle  hits  the  surface,  it  sticks  to  it,  then  a  wedge  having  a 
profile  identical  to  that  of  the  gas  particle  path  might  generate  a  straight 
oblique  shock  wove.  It  should  be  noted,  however,  that  in  the  case  of  a  sticky 
wedge  surface,  solid  particles  are  drawn  away  from  the  suspension  and 
hence  their  loading  ratio  q  could  be  affected. 

The  dependence  of  the  equilibrium  deflection  angle,  0q^,  on  the  flow 

nech  number.  Mo.  is  shown  in  figure  9  for^  =  600,  q  =  i,  D  =  200|im, 

C  =  lOOOJ/(Kgm  K>  and  o  =  l500lCgm/m^.  it  is  evident  from  this  figure  that 
the  equilibrium  deflection  angle  increases  as  the  flow  rtach  number 
increases.  However,  as  the  flow  Mach  number  reaches  high  values,  the 
equilibrium  deflection  angle  is  seen  to  approach  an  upper  limit. 
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It  is  a  common  practice  to  use  shock  polars  to  study  oolique  shock 
waves.  Thus,  in  the  following,  shock  polars  in  dusty  gas  suspension  are 
presented.  • 

The  (P.0g)-polar  for  the  case  Mo  =  2.5,  D  =  lOOpm,  C  =  looo  j/(Kgm 
K),  0  =  1500Kgm/m3  and  various  values  of  q  are  shown  in  figure  to.  The 
(P.Bg)-polar  in  this  case  represents  the  conditions  at  the  end  of  the 
relaxation  zone,  i.e.,  P  =  Peq  and  8g  s  8g^.  As  can  he  seen  the  equilitnum 

pressure,  Peq,  becomes  higher  when  q  is  increased.  In  addition,  the  maximum 
deflection  angle,  known  os  the  detachment  angle,  is  also  seen  to  increase 
with  increasing  q.  While  for  a  dust-free  (pure)  gas  (i.e.,  q  =  0)  >s 

about  29^,  it  reaches  a  value  of  about  56^  when  the  loading  ratio  is  q  =  2. 

A  typical  (^,  8g)-poloris  shown  in  figure  1 1  for  the  some  conditions 
of  figure  tO.  Here  again  Sg  =  Bg  .  It  is  evident  from  this  figure  that  if  a 

certain  deflection  is  required,  say  20^,  then  as  the  loading  ratio  increases 
the  angle  of  incidence,  which  is  needed  to  achieve  the  required  deflection, 
decreases. 

Conclusions 

The  flow  field  developed  when  a  supersonic  dusty  gas  suspension 
passes  through  a  straight  oblique  shock  wave  has  been  investigated. 

The  investigation  included  a  formulation  of  the  governing  equations  of 
the  flow  field  at  hand  and  a  numerical  investigation  of  the  dependence  of 
the  post  shock  suspension  properties  on  the  various  physical  properties  of 
the  solid  particles. 
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Figure  i:  The  shock  wave  sirvjclure  in  a  pure  (solid  lines)  and  a  ousiy 
(dashed  line)  gas. 
a)  lemperalure  profiles 
0)  velocity  profiles 
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Figure  2;  A  schemaUcal  iUustraUon  o(  an  oolique  shock  wove  in  o 
gas  and  the  definition  of  the  (x,y)-  and  (n.s)'Planes. 
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Figure  3:  The  change  of  the  flow  defleciion  angie  with  normal  aisiancc 
from  the  shock  wave  for  various  loading  ratios  q  and  no  =  3, 

^  =  300,  D  =  100  >im,  and  C  =  1000  J/CKgm  K). 
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Figure  4.  The  equilibrium  flow  deflection  angle,  as  a  function  oi  the 
loading  ratio,  q,  for  No  =  4,  ♦  =  600,  D  =  1 00  >im,  end 
C  =  1000  J/tKgmK). 
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Figures:  Tne  change  of  Ihe  flow  deflection  ongle  with  normol  oislonce 
from  the  shock  wave  for  various  solid  panicle  diameters.  D, 

Mq  =  3.  C  =  IOOOJ/(Kgm  K),  ^0  =  300,  and  q  =  0.5. 
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Figure  6:  The  change  of  the  flow  deflection  angle  with  normal  distance 
from  the  shock  wove  for  various  values  of  specific  heat  capacity 
of  the  solid  particles.  C.  end  Ho  s  3.  *0  =  300.  D  =  lOO  )irn.  and  n  = 

0.5. 
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Figure  7;  The  change  of  Ihe  flow  deflection  angle  with  distance  from  me 
shock  wave  for  various  values  of  malerial  density  of  the  solid 
particles,  o,  and  Mo  s  3,  =  300.  0  =  lOO  C  =  1000  J/(Kgm 

K),  and  t\  =  0.5. 
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Figure  6;  The  gas  ana  solid  particles  paths  in  the  (x,y}-plane  for  no  =  3. 
^0  •  60^/  D  s  500  C  :  1000  J/(Kgm  K),  os  1000  Kgm/m^,  and 
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i  FLOW  MACH  NUMBER -Mq 

Figure  9:  The  dependence  of  the  equilibhum  flow  deflection  angle.  6^^^.  on 

the  incident  flow  Mach  number,  for^  s  60<}.  0  =  200 
C  =  1000  J/(Kgm  K).os  1000  Kgm/m^  and  n  =  i. 
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The  Flow  of  a  Dust-Gas  Suspension  with  Dust 
Having  Various  Distributions  of  its  Physical  Properties 
Through  a  Normal  Shock  Wave. 


Abstract 

The  effect  of  distributions  of  non-uniformities  in  the  diameter,  heat 
capacity  and  material  density  of  small  spherical  particles  on  the  properties 
of  a  dust-gas  suspension  passing  through  a  normal  stationary  shock  wave  is 
studied  numerically. 

It  was  found  that  the  gas  temperature  is  practically  independent  of 
the  size  distribution  of  the  dust  particles.  The  suspension  pressure, 
however,  is  very  sensitive  to  the  size  distribution.  It  rises  very  sharply 
when  the  suspension  contains  mainly  small  particles. 

A  distribution  in  the  specific  heat  capacity  of  the  solid  particles 
results  in  a  minor  effect  on  the  velocities  and  the  pressure,  and  a 
pronounced  effect  on  the  temperature  of  the  solid  particles. 

A  distribution  in  the  material  density  of  the  solid  particles  affects 
both  the  velocities  and  the  temperatures  of  the  solid  and  gaseous  phases. 


Introduction 

The  interest  in  the  gas-dynamic  behaviour  of  a  gas-particle 
suspension  grew  in  the  past  three  decades  due  to  its  application  in  many 
engineering  problems.  Some  typical  examples  are;  metalized  propellants  of 
rockets,  jet-type  dust  collectors  and  blast  waves  in  dusty  atmospheres.  In 
addition,  mixtures  with  gases  heavily  laden  with  particles  occur  frequently 
in  industrial  processes  such  as  plastics  manufacturing,  flour  milling,  coal- 
dust  conveying,  powder  metallurgy  and  powdered-food  processing.  General 
descriptions  of  such  flows  can  be  found  in  several  books  and  review  papers 
iSoo  ( 1 967),  Marble  ( 1 970)  and  Rudinger  ( 1 973)1 

The  major  differences  between  the  flow  fields  which  are  developed 
behind  a  normal  shock  wove  in  a  dusty-gos  and  a  pure  (dust-free)  gas  are 
illustrated  in  figures  la  and  b  for  the  temperatures  and  the  velocities. 


respectively.  When  a  steady  pure  gas  encounters  a  normal  shock  wave  it 
experiences  a  sharp  (almost  discontinuous)  change  in  its  thermodynamic  and 
kinematic  properties.  This  sudden  change  is  shown  in  figure  1  to  occur 
between  (0)  and  (1).  The  thickness  of  this  disturbance.  If.  is  only  a  few 
mean  free  paths  of  the  gas  atoms  or  molecules.  Beyond  (1)  the  gas 
properties  remain  constant  (solid  lines  in  figures  la  and  b)  provided  the  gas 
conditions  at  (1)  are  not  sufficient  to  excite  the  internal  degrees  of 
freedom  of  the  gas. 

If,  however,  the  gas  is  laden  with  solid  particles  then  the  suspension 
which  was  originally  at  a  state  of  thermodynamic  and  kinematic 
equilibrium,  ahead  of  the  shock  front,  is  suddenly  changed  into  a  non¬ 
equilibrium  state,  because  the  solid  particles,  due  to  their  size  cornpareo 
with  If,  do  not  experience  any  noticeable  change  in  their  properties  upon 
moving  from  (0)  to  (1).  Thus,  at  (I)  the  gas  has  a  much  higher  temperature 
than  the  dust,  Tg  >  >  Tp^  and  a  much  lower  velocity,  u  <  <  v  (u  is  the  gas 
velocity  and  v  is  the  velocity  of  the  solid  particles).  Morgenthaler  (1962) 
indicated  that  this  is  true  even  if  the  particle  diameter  is  as  small  as  0.1 
)im  (for  shock  waves  in  air  at  nearly  standard  conditions,  where  the  mean 
free  path  is  about  0.066  )im).  Therefore,  the  particles  are  not  influenced  by 
the  initial  disturbance,  and  the  gas  properties  at  (1)  can  be  safely  assumed 
to  be  identical  to  those  of  a  pure  gas  with  the  same  initial  conditions. 

Far  downstream  of  (1),  i.e.,  at  («>)  in  figure  1,  the  gas  and  the  solid 
phases  reach  a  new  state  of  thermodynamic  and  kinematic  equilibrium  via 
momentum  and  energy  exchange.  Theoretically  all  shock  waves  in  dusty 
gases  are  infinitely  thick,  since  equilibrium  is  approached  asymptotically. 
However,  it  is  a  common  practice  to  assign  to  the  shock  wave  an  effective 
thickness  which  is  defined  by  a  requirement  that  the  suspension  properties 
come  close  to  their  equilibrium  downstream  values.  It  was  shown  by 
Gottlieb  &.  Coskunses  (1985)  that  the  suspension  equilibrium  properties  (at 
infinity)  can  be  calculated  from  the  usual  normal  shock  wave  relations, 
provided  that  the  usual  pure  gas  parameters  y  end  R  (the  specific  heat 
capacities  ratio  and  the  specific  gas  constant)  ore  replaced  by  effective 
values  y  and  R  which  solely  depend  on  the  initial  conditions  of  the 
suspension.  This  effective  thickness  is  known  in  the  literature  as  the 
relaxation  zone,  for  it  is  analogous  to  the  relaxation  zone  in  pure  gases 
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where  the  internal  degrees  of  freedom  are  excited.  The  extent  of  the 
relaxation  zone  strongly  depends  on  the  momentum  and  heat  transfer 
mechanisms  which  enable  the  solid  and  the  gaseous  phases  to  reach  a  new 
equilibrium  state.  The  analysis  of  the  relaxation  zone  was  studied  by  many 
investigators.  The  pioneering  works  of  Carrier  (1958),  Kriebel  (1964)  and 
Rudinger  (1964)  verified  the  existence  of  this  relaxation  zone  and  identified 
the  parameters  affecting  it,  namely;  the  solid  particle  diameter,  D,  its  heat 
capacity,  C,  its  material  density,  a,  and  the  loading  ratio,  r\.  Igra  &  Ben-Dor 
(1980)  compared  various  correlations  for  the  drag  coefficient,  Co,  and  the 
heat  transfer  coefficient,  Nu,  and  pointed  out  their  effect  on  the  extent  of 
the  relaxation  zone.  In  addition  they  studied  the  role  of  thermal  radiation 
heat  transfer  between  the  two  phases  and  showed  that  it  can  be  neglected 
when  the  incident  shock  wave  Mach  number  is  smaller  than  five. 

In  all  the  above  mentioned  works,  as  well  as  in  many  others,  the 
gaseous  phase  was  assumed  to  behave  as  a  perfect  gas.  This  assumption 
was  recently  relaxed  by  Ben-Dor  &  Igra  (1982)  and  Igra  &  Ben-Dor  (1984) 
who  solved  the  flow  field  while  accounting  for  real  gas  effects. 
Dissociating  nitrogen  was  the  gaseous  phase  in  the  latter  work  and  ionizing 
argon  in  the  former. 

The  assumption  that  the  solid  particles  are  inert,  which  was  also 
adopted  in  most  of  the  published  studies,  was  also  relaxed  by  Elperin,  Ben- 
Dor  &  Igra  (1986)  who  solved  the  flow  field  of  an  oxygen-carbon  suspension 
passing  through  a  normal  shock  wave,  behind  which  the  carbon  particles 
reached  their  ignition  temperature  and  burned  out. 

In  the  present  study,  another  commonly  used  assumption  concerning 
the  solid  phase  is  relaxed.  It  has  been  a  common  practice  to  assume  that 
the  solid  phase  consists  of  spherical  particles  of  identical  diameter,  D, 
identical  material  density,  a,  and  heat  capacity,  C.  In  reality,  however,  this 
is  not  the  case,  for  the  particles  comprising  the  solid  phase  do  not  have,  in 
general,  a  uniform  size  even  though  they  might  have  the  same  physical 
properties,  i.e.,  a  and  C.  Furthermore,  there  could  be  cases  in  which  the 
suspension  consists  of  different  solid  materials,  i.e.,  a  and  C. 

The  purpose  of  this  study  is  to  numerically  investigate  the  influence 
of  the  above  mentioned  solid  phase  non-uniformities  on  the  flow  field  which 
is  developed  behind  a  normal  shock  wave. 
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In  the  following,  the  basic  assumptions  upon  which  the  present  model 

is  based  are  given.  The  assumptions  are  followed  by  the  governing  equations 

and  the  numerical  results  arising  from  their  solution. 

Theoretical  Background 

% 

Assumptions 

The  assumptions  upon  which  the  present  model  is  based  and  their 

implications  are  listed  in  the  following: 

1)  The  gaseous  phase  behaves  as  an  ideal  gas.  Thus,  the  equation  of  slate 
of  the  gas  is  P  =  pgRTg.  Note  that  it  is  not  assumed  here  that  the  gas  is 
calorically  ideal.  Alternatively,  the  dependence  of  both  Cp  and  Cy  on 
the  gas  temperature  is  accounted  for.  This  has  not  been  done  in 
previous  studies  where  both  Cp  and  Cy  were  assumed  to  be  constant. 

2)  All  solid  particles  are  rigid,  chemically  inert  small  spheres  uniformly 
distributed  in  the  gaseous  phase.  Thus  there  is  no  heat  addition  or 
reduction  due  to  chemical  processes  between  the  solid  and  the  gaseous 
phases.  Furthermore,  Re  and  Nu  ore  based  on  the  particle  diameter,  D. 

3)  The  volume  of  the  solid  phase  in  the  suspension  can  be  neglected.  Thus 
the  momentum  and  energy  exchange  between  the  solid  particles  can  be 
ignored. 

4)  Aside  from  momentum  and  energy  interactions  between  the  gaseous  and 
the  solid  phases,  the  gaseous  phase  is  considered  to  be  a  perfect  fluid, 
i.e.,  no  other  viscous  or  conduction  effects  are  considered. 

5)  The  solid  particles  are  too  large  to  experience  any  change  in  their 
thermodynamic  and  dynamic  properties  upon  their  passage  through  the 
shock  front.  In  addition  they  are  also  large  enough  not  to  experience 
Brownian  motion.  Thus,  the  partial  pressure  of  the  solid  phase  can  be 
neglected. 

6)  The  solid  particles  are  small  enough  to  satisfy  the  condition  B,  <0  1. 
where  Bj  is  the  Biot  number.  Bj  =  hr/kp  (h  is  the  coefficient  of  heat 
transfer,  r  is  the  radius  of  the  particle,  and  kp  is  its  thermal 
conductivity).  Thus  the  temperature  within  the  solid  particles  can  be 
assumed  to  be  uniform. 
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7)  The  weight  of  the  solid  particles  and  the  buoyancy  forces  experienced 
by  them  are  negligibly  small  in  comparison  with  the  drag  forces  acting 
on  them. 

8)  The  heat  capacity,  C,  of  the  solid  particles  is  constant. 

9)  Ahead  of  the  normal  shock  wave  the  suspension  is  at  a  state  of 
thermodynamic  and 'dynamic  equilibrium,  i.e.,  uq  =  vq  and  Tgg  =  Tpg, 

where  u  and  v  are  the  velocities  of  the  gas  and  the  solid  particles,  and 
Tg  and  Tp  are  the  temperatures  of  the  gas  and  solid  particles, 
respectively. 

in  addition  to  the  above  listed  assumptions  it  is  assumed  that  the  flow  field 
under  consideration  is  one-dimensional  and  steady. 

Governing  Equations 

Prior  to  listing  the  governing  equations  it  should  be  mentioned  that 
the  considered  suspension  is  composed  of  n  +  1  phases.  One  phase  is  the 
gaseous  phase.  The  remaining  n  phases  are  all  solid  phases.  Each  of  the  n 
solid  phases  consists  of  identical  particles.  However,  each  one  of  the  n 
solid  phases  differ  from  all  the  other  solid  phases  in,  at  least,  one  of  the 
following  properties;  the  diameter  of  the  solid  particle,  D,  its  material 
density,  a,  or  its  heat  capacity,  C. 

Based  on  the  foregoing  assumptions  the  governing  equations 
describing  the  problem  at  hand  are: 

-  continuity  of  the  gaseous  phase; 

-^(pqu)  =  0  (1) 

-  continuity  of  the  i-th  solid  phase; 


^(ppjVj)  =  0  i=lton 


(2) 
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-  conservation  of  momentum  of  the  gaseous  phase; 

-|^(p,u2»P)  =  |FDj  (3) 

-  conservation  of  momentum  of  the  i-th  solid  phase; 

(ppi  )  = -Fpi  (4) 

-  conservation  of  energy  of  the  gaseous  phase; 

PgU  -|^[CpTg  ^  1/2  u2l  =  Vi  (5) 

-  conservation  of  energy  of  the  i-th  solid  phase; 

Ppi  Vi  iCi  Tpi  +  1/2  v^l  =  -Qpi  -  Fpi  Vi  (6) 

-  equation  of  state  of  the  gaseous  phase; 

P  =  pqRTg  (7) 

In  the  above  equations  pg,  u,  Tg  and  P  are  the  density,  velocity,  temperature 
and  pressure  of  the  gaseous  phase,  respectively,  ppj,  Vi  and  Tpi  are  the 
spatial  density,  velocity  and  temperature  of  the  I-th  solid  phase, 
respectively.  Note  that  the  material  density  of  the  solid  particles  is  oj 
which  can  be  related  to  the  spatial  density  ppi  via  the  relation  ppj  =  ni  Vpj  oj 
where  ni  is  the  number  density  of  the  i-th  solid  particle  and  Vpj  is  the 
volume  of  a  single  i-th  solid  particle.  It  should  also  be  noted  that  according 
to  assumption  5  the  partial  pressures  of  all  the  solid  phases  are  zero.  For 
this  reason  P  is  not  only  the  pressure  of  the  the  gaseous  phase,  but  the 
pressure  of  the  suspension  as  well.  Cp  and  R  are  the  specific  heat  capacity 
at  constant  pressure  and  the  specific  gas  constant  of  the  gaseous  phase  ana 
Cj  is  the  heat  capacity  of  the  particles  of  the  i-th  solid  phase.  Fq,  is  the 
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drag  force  per  unit  volume  exerted  by  the  gaseous  phase  on  the  particles  of 
the  i-th  solid  phase  and  Qpj  is  the  heat  transferred  per  unit  volume  from  the 
particles  of  the  i-th  solid  phase  to  the  gaseous  phase. 

The  drag  force  Fpi  can  be  calculated  from: 

Fqj  =3/4  pgPpj  (vj  “  u)  I  Vi  -  u  I  Cp^  /  (Dj  ff])  (3) 

where  Di  is  the  diameter  of  the  solid  particle,  a\  is  their  material  density, 
and  CD|  the  drag  coefficient  is  given  by 

COj  =  0.48  +  28  Rer°  ®5  (9) 

The  slip  Re  number  for  the  i-lh  solid  phase  is  given  by 


Rej 


Pgivi  -  u  ID] 
Fg 


(10) 


where  the  dynamic  viscosity,  pg,  depends  on  the  gas  temperature,  Tg. 
The  transferred  heat,  Qpj,  can  be  calculated  from: 


Qpj  =  6hj  ppj  (Tpj  -  Tg)  /  (D^  (jj) 


(11) 


where  hj,  the  coefficient  of  convection  heat  transfer  can  be  obtained  from 


Nui  is  the  Nusselt  number  of  the  i-th  solid  phase.  It  can  be  obtained  from 

NUi  =  2 0.459  PrO-33  RejO-SS  (13) 


In  the  last  expression  Pr  is  the  Prandtl  number; 


Pr  = 


Me 

kg 


(14) 


where  kg  (also  in  12)  is  the  thermal  conductivity  of  the  gaseous  phase. 
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The  above  set  of  governji.g  equations  consists  of  4  +  3n  equations 
(conservation  of  mass,  momentum  and  energy  for  the  gaseous  phase, 
equation  of  state  of  the  gaseous  phase  and  conservation  of  mass,  momentum 
and  energy  for  each  of  the  n  solid  particles).  The  number  of  the  unknowns  in 
the  above  set  of  equations  is  also  4+3n  (pg,  u.  Tg  and  P  for  the  gaseous  phase 
and  ppj,  Vj  and-Tpj  for  each  of  the  n  solid  phases).  Thus  the  above  set  of 
governing  equations  is  solvable  in  principle. 

Numerical  Approach 

Computer  Code 

There  are  many  computational  packages  capable  of  numerically 
solving  differential  equations.  The  one  chosen  for  solving  the  differential 
equations  governing  the  problem  at  hand  is  the  ItlSL  (International 
Mathematical  &  Statistical  Libraries)  package  suited  to  operate  on  a  CDC 
Cyber  840  computer.  The  IMSL  package  contains  three  computer  codes  for 
solving  differential  equations  with  given  initial  conditions;  they  are; 

DVERK  -  based  on  the  Runge-Kutta  method,  and  recommended  for  cases 
where  high  accuracy  is  not  required  and  where  the  derivatives  can 
be  simply  calculated. 

DGEAR  -  based  on  the  predictor-corrector  method.  Although  it  results  in 
poor  accuracy  it  is  preferable  to  DVERK  when  the  calculation  of 
the  derivatives  is  difficult  and  hence  expensive. 

DREBS  -  based  on  the  extrapolation  method.  Preferable  when  high 
accuracies  are  required  and  when  the  derivatives  can  be 
calculated  in  a  relatively  simple  and  inexpensive  way. 

The  DREBS  code  was  adopted  in  the  present  study  because  it  is  the 
most  accurate  one.  This  code  performs  a  triple  check  of  the  obtained  error 
between  each  two  extrapolation  steps.  Since  the  DREBS  code  is  limited  to 
the  case  where  each  derivative  is  independent  of  the  other  derivatives  the 
governing  equations  need  to  be  rearranged.  The  rearrangement  of  the 
governing  equations  (1)  to  (7)  results  in; 


dx 


FpiVj  - 


Fpi  u3^ 


fr 

■  u2-RTg^ 


(15) 
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du  ~  Pi'* 

dx  “  u-RTfl 

^  pf-^Ea.  To  du) 
dx  =  P5'^[  dx  '  u  dx 

dpq  pq  du 


dx 

u  dx 

dTpi 

Qoi 

dx 

Cj  ppi  Vj 

dvj 

Foi 

dx  " 

Ppi  Vj 

dppi 

dx 

ppi  dVj 

Vi  dx 

(16) 

(17) 

(18) 

(19) 

(20) 

(21) 


Note  that  although  equations  (15)  to  (21)  do  not  reflect  the  dependence  of 
the  gas  specific  heat  capacity  on  its  temperature  since  a  space  gradient  of 
Cp  does  not  appear,  the  value  of  Cp  was  calculated  at  any  position  x 
according  to  the  temperature  of  that  position.  This  procedure  assumes  that 
the  space  gradient  of  Cp  is  neglectable. 

Initial  Conditions 

As  mentioned  earlier  the  properties  of  the  solid  phases  do  not  change 
when  they  pass  through  the  normal  shock  wave.  Thus,  the  solid  particles 
velocity,  v^,  temperature  Tpj,  and  spatial  density,  ppj,  immediately  behind 
the  shock  front  remain  identical  to  their  appropriate  values  ahead  of  the 
shock  front,  i.e., 


(Vi)l  =  (Vi)Q 

(Tpj))  =  (Tpj)o 
(ppj)l  =  (ppi)o 

Furthermore,  the  multi-phase  suspension  is  assumed  to  be  in 
thermodynamic  and  kinematic  equilibrium  ahead  of  the  shock  front. 
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therefore;  (vi)o  =  uo  and  (Tpi)Q  =  To  and  (ppi)o  =  rjjpgg  where  tw  is  the 

loading  ratio  of  the  i-th  solid  phase. 

The  properties  of  the  gas,  on  the  other  hand,  change  almost 
discontinuously  across  the  shock  wave.  The  gas  properties,  immediately 
behind  the  shock  front,  can  simply  be  calculated  from  the  normal  shock  jump 
conditions,  i.e.; 


P]  ^^^0  Y~^ 

Pq  '  '(*■  1  y+ 1 


(23)? 


Pg,  yiWp 
Pgo’ 


'•Jl 


Y"  1  2V  2  ' 


T, 


50 


n 


n,= 


m2  . 

n  + 

_|o _ ^ 

2y  ^2  , 
- T 

y- 1  0 


(y>1)2 
0  2(y-1) 

2 


01  =  Ml  V  Y  R  Tgj 


(24) 


(25) 


(26) 


(27) 


Mo  and  M]  in  the  above  equations  are  the  flow  Mach  numbers  immediately 
ahead  and  behind  the  shock  front. 


The  Final  Conditions 

As  mentioned  the  extent  of  the  relaxation  zone  in  the  case  of  a  dusty 
gas  is  infinitely  long  for  the  equilibrium  values  are  approached 
asymptotically.  However,  as  shown  by  Gottlieb  5^  Coskunses  (1965)  the 
suspension  equilibrium  properties  at  infinity  can  be  calculated  very  simply 
by  using  the  above  listed  normal  shock  wave  equations  [(23)  to  (27)1  and 


replacing  Y  by  y  and  R  by  R.  In  this  approach  the  multi-phase  suspension  is 
assumed  to  be  a  single  phase  gas  having  new  values  for  y  and  R.  Based  on 

Gottlieb  &.  Coskunses  (1985),  the  relations  for  y  and  R  are  as  follows; 


-  ^Hi^i  •^  (l~Sri,)Cp 

+  (l-2Tij)Cv 


(28) 


R  =  (l-2r|j)R 


(29) 


where 


Ppi 


Ippi  +  Pg 


(30) 


Note  that  in  many  papers  the  loading  ratio  is  defined  as  t^.  =  pp/pg.  Thus, 
the  two  different  definitions  can  be  related  by  the  following  expression: 


Ti  = 


2ri.  +  1 


(31) 


(It  is  clear  from  these  relations  that  if  the  total  loading  ratio  Zq,  and  heat 

capacity  C\  are  constant,  then  y  and  R  are  constant  and  hence  the 
equilibrium  properties  at  the  end  of  the  relaxation  zone  are  identical.)  These 
new  values,  i.e.,  y  and  R  yield  a  new  speed  of  sound  a.  which  satisfies; 


ao  <  aq  =  (yRTgQ)i'^2 


Thus  the  flow  Mach  number  ahead  of  the  shock  wave  becomes  no  >  no.  By 
replacing  y,  R  and  ho  in  equations  (23)  to  (27)  with  y  ,  R  and  ho,  one  can 
calculate  the  equilibrium  properties  of  the  suspension  at  infinity. 

The  knowledge  of  the  suspension  properties  at  infinity,  i.e.,  the 
equilibrium  final  conditions  can  serve  as  an  excellent  means  of  checking  the 
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reliability  of  the  numerical  results  since  they  should  be  approached 
asymptotically. 

Since  the  equilibrium  conditions  at  infinity  are  known  apriori  the 
integration  was  terminated  when  the  suspension  properties  came  as  close 
as  2%  of  their  corresponding  equilibrium  value.  The  distance  where  the 
velocities  reached  this  condition  is  called  the  kinematic  relaxation  length, 
ly,  and  the  distance  where  the  temperatures  reached  this  condition  is  called 
the  thermal  relaxation  length,  Ij.  These  two  relaxation  lengths  are  in 
general  different. 

Phusical  Properties  of  the  Various  Phases 

The  gaseous  phase  was  chosen  to  be  nitrogen.  Yum,  Weissman  &  Mason 
(1962)  have  calculated  the  transport  properties  of  nitrogen  and  presented 
their  results  in  a  tabular  form  for  both  the  dynamic  viscosity,  and  the 
thermal  conductivity,  kg.  A  least  square  polynomial  fit  to  their  results 
yields  the  following  expressions  for  kg  andjig; 

kg[W/(mK)l  =  1 .386  xl 0-2  5.3 1 1  x  1 0'S  Tg  -  9.822  x  1 0*  1 0  T^ 
and 

MglKgmArn  see))  =  0.7226  x  10-5  +  2.768  X  lO'Q  Tg-0.5933  x  10->2t^ 

The  expression  for  the  specific  heat  capacity  at  constant  pressure,  Cp,  was 
adopted  from  Sontag  &  Van-Wylen  (1971); 

Cnlm2/(sec2  K)1  =  1.3686X  103-  2.786x  IQST’’  +5.1741  x  10?T’’ 

H  g  g 

In  the  above  expression  Tg  is  in  degree  K. 

Five  solid  phases  were  chosen  in  order  to  numerically  simulate  the 
effect  of  non-uniformities  of  the  solid  particles  on  the  flow  field.  Thus,  in 
equations  (18)  and  (20)  n  =  5  and  each  of  equations  (15)  to  (17)  represents 
five  equations  with  i  =  1,  2, ...,  5. 

Each  solid  phase  is  identified  by  three  parameters.  The  diameter  of 
its  particles,  Dj,  the  heat  capacity  of  its  particles  Cj  and  their  material 
density,  oj. 
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Nine  different  cases  (for  Mo  =  1.5)  were  investigated.  They  are  shown 
in  tables  1  to  3.  In  table  1  all  the  five  solid  phases  have  the  same  material 
density  ci  =  az  =•••  =  as  =  a  and  heat  capacity  Ct  =  C2  =  -  =  C5  =  C  but 
different  particle  diameters  Di  =  5jim,  D2  =  7.5pm,  D3  =  10pm,  D4  =  25pm 
and  Ds  =  50pm.  The  seven  cases  differ  in  the  loading  ratios  of  the  various 
solid  phases. 

The  partial  loading  ratios  r|i  are  chosen  in  such  a  way  that  the  total 

5 

loading  ratios  for  each  case,  i.e,  qtoW  =  Xni  eQuals  always  0.2. 

1 

Thus,  for  all  the  seven  cases,  the  suspension  properties  should  reach 
the  same  equilibrium  values  at  the  end  of  the  relaxation  zone. 

In  the  case  shown  in  table  2  also  for  Mo  =  1.5,  the  five  solid  phases 
have  identical  particle  diameters,  D  =  50pm  and  material  densities,  a  = 
1500  Kgm/m3,  but  different  heat  capacities  Ci  =  100,  Cz  =  500,  C3  =  1000, 
C4  =  5,000  and  C5  =  10000  J/(Kgm  K).  The  loading  ratio  of  each  solid  phase 
is  q  =  0.04  and  therefore  the  total  loading  ratio  is  again  0.2. 

The  last  case  is  shown  in  table  3.  Here  all  the  solid  phases  have 
identical  particle  diameters,  D  =  lOpm  and  heat  capacities,  C  =  1000 
J/(Kgm  K),  but  different  material  densities,  a\  -  1000,  a2  =  1250,  03  = 
1500,  04  =  1750  and  as  =  2000  Kgm/m^.  The  incident  shock  wave  Mach 
number  is  again  Mq  =  1.5  and  the  total  loading  ratio  is  0.2  (qj  =  0.04  for  each 
solid  phase). 

Numerical  Results  and  Discussions 
Distribution  of  the  Diameter  of  the  Dust  Particles 

Figures  2  to  7  illustrate  the  flow  field  properties  which  are 
developed  behind  a  normal  shock  wave  with  Mq  =  1.5  in  cases  1  to  6, 
respectively  (see  table  1).  Each  of  these  six  figures  contains  six  velocity 
and  six  temperature  profiles  (one  for  the  gaseous  phase  labelled  with  "G" 
and  five  for  each  of  the  five  solid  phases  labelled  "1"  to  "5“  in  accordance 
with  table  1)  and  one  pressure  profile.  The  dashed  lines  in  the  right  hand 
sides  of  figures  2  to  7  indicate  the  equilibrium  values  at  infiniiy,  which 
should  be  approached  asymptotically.  As  mentioned  earlier  these  values  are 
calculated  apriori. 


It  is  clear  from  figures  2  to  7  that  the  larger  the  particle  diameter 
is,  the  longer  it  takes  for  it  to  reach  a  dynamic  equilibrium.  This  is  due  to 
the  fact  that  the  inertial  forces  depend  on  while  the  drag  forces  are 

2 

roughly  proportional  to  .  Thus  the  effectiveness  of  the  drag  forces  in 

slowing  down  the  solid  particles  decreases  as  the  diameter  of  the  solid 
particles  increases.  The  solid  phase  with  the  larger  particles  is  also  the 
last  to  reach  thermal  equilibrium  with  the  gaseous  phase.  This  is  not 
surprising  since  the  larger  the  solid  particle  is,  the  slower  will  be  its 
temperature  rise.  Thus  the  solid  phase  with  the  smallest  particles  (D  = 
5jim)  reaches  a  thermal  equilibrium  with  the  gaseous  phase  within  few 
centimeters  while  it  takes  tens  of  centimeters  for  the  solid  phase  with  the 
largest  particles. 

The  velocity,  temperature  and  pressure  profiles  of  the  gaseous  phase 
for  cases  1  to  6  are  shown  in  figure  8.  It  is  clearly  seen  in  figure  8  that  all 
the  six  cases  approach  asymptotically,  identical  equilibrium  properties  at 
infinity.  This  fact  should  not  be  surprising  for  the  total  loading  ratio  in  all 
the  six  cases  is  g  =  0.2  and  the  heat  capacity  of  the  solid  particles  is  the 
same.  Although  the  equilibrium  properties  are  identical,  the  flow 
properties  inside  the  relaxation  zone  do  depend  on  the  size  distribution  of 
the  particles  of  each  solid  phase.  The  dependence  is  minimal  for  the  gas 
temperature  where  the  six  profiles  merge  into  an  almost  one  curve,  and  it  is 
maximal  for  the  suspension  pressure.  The  suspension  pressure  increases 
very  sharply  when  it  consists  mainly  of  small  particles  (case  5)  and  it 
increases  very  slowly  when  it  consists  mainly  of  large  particles  (case  6). 
Thus  it  can  be  concluded  that  increasing  the  mass  fraction  of  the  small 
solid  particles  in  the  suspension  results  in  a  faster  rise  in  the  suspension 
pressure.  Similarly,  the  suspension  with  the  largest  amount  of  small 
particles  results  in  the  sharpest  decrease  in  the  gaseous  phase  velocity.  The 
suspension  with  the  largest  amount  of  large  particles,  experiences  the 
smallest  decrease  in  the  gaseous  phase  velocity.  It  is  interesting  to  note 
that  the  gaseous  phase  velocity  falls  from  its  value  immediately  behind  the 
shock  front  in  spite  of  the  fact  that  it  exchanges  momentum  with  the  solid 
phases  which  have  much  higher  post  shock  velocities.  The  reason  for  this 
peculiar  behaviour  lies  in  the  fact  that  the  dust  presence  causes  a  very 
large  increase  in  the  density  of  the  gaseous  phase.  This  in  turn  results  in  a 


large  decrease  in  the  gaseous  phase  velocity  since  the  flow  field  is  one 
dimensional  and  p^u  must  remain  constant. 

The  general  shapes  of  the  various  property  profiles,  for  the  six  cases 
of  multi-solid-phases,  are  similar  to  those  obtained  when'  the  suspension 
contains  a  single-solid-phase  only,  i.e.,  case  7  in  table  1,  whose  solution  is 
seen  in  figure  9.  Since  for  this  case  =  0.2  too,  the  equilibrium  properties 
at  infinity  are  identical  to  those  shown  in  figures  2  to  7. 

The  thermal  and  kinematic  relaxation  lengths  for  cases  1  to  7  are 
given  in  table  4.  The  earlier  remark  that  the  temperature  of  the  gaseous 
phase  is  almost  unaffected  by  the  size  distribution  of  the  particles  of  the 
solid  phases  (see  figure  8,  where  all  the  temperature  profiles  merge  into  a 
single  line  at  a  very  short  distance  behind  the  shock  fronts)  is  further 
supported  by  the  fact  that  the  thermal  relaxation  length  It  is  almost  the 
same  for  the  six  multi-phase  cases. 

The  difference  between  the  kinematic  relaxation  lengths,  ly,  is  much 
larger.  The  two  extreme  cases,  5  &  6,  differ  by  about  8.5^.  Note  that  the 
difference  in  the  thermal  relaxation  length  between  the  two  extreme  cases 
is  only  about  1.5^. 

However,  the  difference  in  both  the  thermal  and  kinematic  relaxation 
lengths  between  the  multi-solid  phases  (case  1  to  6)  and  the  single-solid 
phase  (case  7)  is  enormous.  Both  lx  and  ly  are  more  than  90%  smaller  when 
the  suspension  contains  a  single  solid  phase  with  D  =  lOpm.  It  is  obvious 
that  the  enormous  increase  in  both  lx  and  ly  when  the  suspension  contains  a 
size  distributed  solid  phase  arises  from  the  presence  of  the  large  particles, 
D  =  50)im,  in  the  six  cases.  This  is  due  to  the  fact  that  the  large  particles 
slow  down  and  heat  up  very  slowly  in  comparison  with  smaller  particles. 
Thus  it  can  be  concluded  that  the  extent  of  the  relaxation  zone  is 
determined  solely  by  the  solid  particle  having  the  largest  diameter. 

Distribution  of  the  Heat  Capacity  of  the  Dust  Particles 

Case  8  (table  2)  in  which  the  solid  particles  of  the  five  phases  have 
identical  diameters  D  =  50)im  and  material  density  o  =  1500  Kgm/rn^  but 
different  heat  capacities,  C,  were  solved.  The  results  are  shown  in  figure 
10  for  Mo  =  1.5. 

It  is  evident  from  figure  10  that  in  spite  of  the  fact  that  the  heat 
capacity,  C,  of  each  of  the  solid  phases  is  different,  they  all  have  the  same 


velocity  and  pressure  inside  the  relaxation  zone.  Thus,  the  heat  capacity 
has  no  affect  on  the  velocities  of  the  solid  particles  and  the  gaseous  phases 
and  does  not  affect  the  suspension  pressure.  However,  the  five  solid  phases 
have  different  temperature  profiles.  The  temperature  of  the  solid  particles 
with  the  highest  value  of  C  experiences  the  slowest  increase.  For  the 
smallest  value, of  C  the  temperature  rise  of  the  solid  particles  is  extremely 
sharp.  It  overshoots  the  temperature  of  the  gaseous  phase  before  it  drops 
down  to  the  equilibrium  value  at  infinity. 

The  thermal  and  kinematic  relaxation  lengths,  Ij  and  ly,  are  also 
shown  in  table  4  for  this  case  (case  8).  The  large  values  for  Ij  arise  from 
the  fact  that  it  takes  a  long  time  for  the  solid  phase  with  the  largest  heat 
capacity  to  reach  equilibrium  with  the  gaseous  phase  because  their 
temperature  rise  is  very  slow.  The  large  values  for  the  kinematic 
relaxation  length  are  due  to  the  fact  that  the  diameter  of  all  the  solid 
particles  was  chosen  to  be  D  =  50)im  for  this  case. 

Distribution  of  the  Material  Densitu  of  the  Dust  Particles 

Case  9  (table  3)  in  which  the  solid  particles  of  the  five  phases  have 
identical  diameters  0  =  10}im  and  heat  capacities  C  =  1000  J^'Kgm  K)  but 
different  material  densities  is  shown  in  figure  M  for  Mo  =  1.5.  Case  7 
which  has  identical  values  of  D  and  C,  as  well  as  the  total  loading  ratios  of 
q  =  0.2,  but  a  constant  value  of  a  =  1500  Kgm/m^  is  added  to  figure  1 1  in 
dashed  lines. 

The  solid  lines  in  figure  1 1  indicate  that  the  distribution  in  the  solid 
particles  density  results  in  a  decrease  in  the  gas  pressure,  temperature  and 
velocity  gradients  behind  the  shock  front.  The  equilibrium  values,  however, 
at  infinity  are  identical  for  the  two  cases.  In  addition  it  is  evident  from 
figure  1 1  that  the  larger  the  material  density  of  the  dust  particles  is,  the 
slower  its  decay  becomes.  This  is  so  because  its  inertia  depends  linearly  on 
the  material  density.  Similarly,  the  larger  the  materiel  density  of  the  dust 
is,  the  slower  its  temperature  rise  becomes,  because  it  has  a  higher  heat 
capacity.  The  effect  of  the  distribution  in  the  material  density  of  the  dust 
particles  on  the  thermal  and  kinematic  relaxation  lengths,  Ij  and  ly, 
respectively,  can  be  seen  by  comparing  cases  7  &  9  in  table  4.  For  the 
tio  =  1.5  case  (7  &  9)  both  It  and  ly  increase  by  about  25^. 
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Conclusions 

The  flow  field  which  is  developed  behind  a  steady  normal  shock  wave 
in  a  dust-gas  suspension  has  been  solved  numerically  while  accounting  for 
distributions  of  non-uniformities  in  the  physical  properties  of  the  solid 
phase.  The  model  was  bared  on  a  multi-phase  system  which  consists  of  one 
gaseous  phase  and  five  solid  phases.  The  five  solid  phases  differed  from 
each  other  by  one  physical  property  only,  which  was  either  the  diameter  of 
the  solid  particles  or  their  heat  capacity  or  their  material  density. 

For  a  fixed  total  loading  ratio  it  was  found  that  the  gas  temperature 
profile  downstream  of  the  shock  wave  is  practically  independent  of  the 
size  distribution  of  the  dust  particles.  This  might  also  be  a  result  of  the 
relatively  small  loading  ratio  which  was  used  in  the  calculation. 
Consequently,  it  should  not  necessarily  hold  for  an  arbitrary  loading  ratio. 
The  suspension  pressure  profile,  however,  is  very  sensitive  to  the  size 
distribution.  The  suspension  pressure  rises  very  sharply  when  the 
suspension  contains  mostly  small  particles.  Similarly  the  suspension  with 
the  largest  amount  of  small  particles  results  in  the  sharpest  decrease  in 
the  velocity  of  the  gaseous  phase. 

A  distribution  in  the  heat  capacity  of  the  solid  particles  resulted  in  a 
minor  effect  on  the  velocity  of  the  solid  particles  and  practically  no  effect 
on  the  suspension  pressure.  The  effect  on  the  temperature  of  solid  particles 
was  much  more  pronounced.  As  expected,  the  particles  with  the  smaller  heat 
capacity  were  heated  faster  and  their  temperature  reached  that  of  the 
gaseous  phase  earlier  than  those  having  larger  heat  capacities. 

A  distribution  in  the  material  density  of  the  solid  particles  resulted 
in  different  velocity  and  temperature  profiles  for  each  of  the  five  solid 
phases,  and  a  slower  rise  in^he  suspension  pressure. 

Acknowledgement 

The  financial  support  received  from  the  L).  S.  Air  Force  under  Grant 
AFOSR-86-0349  is  acknowledged  with  thanks. 

List  of  References 

Ben-Dor,  G.  &  Igra.  0.,  1982,  J.  Plasma  Phys.,  27,  377. 

Carrier,  G.  F.,  1950,  J.  Fluid  Mech.,  4,  376. 

Elperin,  I.,  Ben-Dor,  G.  &  Igra,  0.,  1986,  ASME  J.  Fluid  Eng.,  108,  354. 


18 


Gottlieb,  JJ.  &  Coskunses,  C.E.,  1 985,  UTIAS  Rep.  295. 

Igra,  0.  &  Ben-Dor,  6.,  1980,  IsraelJ.  Tech.,  18.  159. 

Igra,  0.  &  Ben-Dor,  6.,  1 984,  J.  Plasma  Phys.,  31,  115. 

Kriebel,  A.R.,  1 964,  J.  Basic  Eng.,  Trans.  ASHE,  D85,  655 
Marble,  F.E.,  1970,  Ann.  Rev.  Fluid  Mech.,  2,  397. 

Morgenthaler,  J.H.,  1962.  in  Progress  in  Astronautics  and  Rocketry  Ed.  S.  S. 

Penner&  F.  A.  Williams,  Vol.  6,  Academic  Press,  New  York,  N.Y.,  p.  145. 
Rudinger,  G.,  1964,  Phys.  Fluids,  7,  658. 

Rudinger,  G.,  1973,  Appl.  Mech.  Rev.,  26,  273. 

Sontag,  R.E.  &  Van  Wylen,  G.  J.,  1971.  Introduction  to  Thermodunamics,  John 
Willey  and  Sons 

Soo,  S.L.,  1967,  Fluid  Dynamics  of  Multiphase  Systems.  Blaisdell  Publishing 
Co.,  Waltham,  Mass.,  U.S.A. 

Yum,  K.S.,  Weissman,  S.  &  Mason,  E.A.,  1962,  Phys.  Fluids,  5,  672. 


19 


Figure  Captions 

Figure  1;  Shock  wave  structure  in  pure  and  dusty  gases 

Figure  2:  Flow  field  for  case  1  of  table  1 

Figure  3:  Flow  field  for  case  2  of  table  1 

Figure  4;  Flow  field  for  case  3  of  table  1 

Figure  5:  Flow  field  for  case  4  of  table  1 

Figure  6:  Flow  field  for  case  5  of  table  1 

Figure  7;  Flow  field  for  case  6  of  table  1 

Figure  8;  Gas  properties  for  cases  1  to  6 

Figure  9:  Flow  field  for  case  7  of  table  1 

Figure  10;  Flow  field  for  case  8  of  table  2 

Figure  1 1;  Flow  field  for  case  9  of  table  3 
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Table  1;  The  loading  ratios,  of  the  78110113  size  particles  for 
a  =  1500  Kgm/m^  and  C  =  1000  J/(KgmK) 
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Table  2;  The  specific  heat  capacities,  C,  of  the  five  various  solid  phases. 

For  all  the  solid  phases 

Mo  =  1 .5,  D  =  50  X  1 0-6m,  (J  =  1 500  Kgm/m3^  n  =  0.04 


Ci(J/(Kgm  K)1 

Case  Number 
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Table  3;  The  dust  material  densities,  oj,  of  the  five  various  solid  phases. 

For  all  the  solid  phases 

Mo=  1.5,  D=  lOx  10-6m,  C=  1000J/(Kgm  K),  n  =  0.04 


ffilKgm/m^l 
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